cosint(5.2) 
cosint([0:0.1:1])
syms x;
f = cosint(x);
diff(x)
